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Abstract 

We consider the formation of a ring-like magnetic field in collisions of bubbles of 
broken phase in an abelian Higgs model. Particular attention is paid on multiple 
collisions. The small collision velocity limit, appropriate to the electroweak phase 
transition, is discussed. We argue that after the completion of the electroweak 
phase transition, when averaged over nucleation center distances, there exists a 
mean magnetic field B ~ 2.0 x 10^'^ G with a coherence length 9.1 x 10^ GeV"^ 
(for niH = 68 GeV). Because of the ring-like nature of B, the volume average 
behaves as i? ~ 1/L. Taking into account the turbulent enhancement of the 
field by inverse cascade, we estimate that colliding electroweak bubbles would 
give rise to a mean field Brms — 10~^^ G at 10 Mpc comoving scale today. 
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1 Introduction 



Cosmological first order phase transitions may give rise to primordial magnetic fields, 
which then could act as the seed field required for the dynamo explanation of the 
observed galactic magnetic fields [0. A first order phase transition proceeds by nucle- 
ation of bubbles of broken phase in the background of unbroken phase. In the case of a 
first order electroweak (EW) phase transition, the Higgs field inside a given bubble has 
an arbitrary phase. The bubbles expand and eventually collide, while new bubbles are 
continuously formed, until the phase transition is completed. This also involves the 
equilibration of the phases of the complex Higgs fields, the gradients of which act as a 
source for gauge fields, thus making the generations of magnetic fields possible. The 
growth of the bubble is much affected by the non-linear interplay between the field 
configuration constituting the bubble and the background plasma. In the EW case 
hydrodynamical studies show ^ that the expanding bubble is preceded by a shock 
front, which heats up the universe back to the critical temperature. The transition is 
then completed by the merging of the slowly expanding bubbles. 

There are two different, but not mutually exclusive, theoretical scenarios for gen- 
erating magnetic fields in first order phase transitions. One employs directly gradients 
of the complex Higgs field in collisions between bubbles of broken phase, as discussed 
in 0. The other is based on the spontaneous appearance of electrical currents and 
turbulent fiow near the bubble walls and has been applied both to QCD and EW 
[^, ^ phase transitions. In the EW phase transition separation of electric charges oc- 
curs at the phase boundary because of baryon number gradients. These give rise to a 
net current and hence magnetic fields, the fate of which is dependent on the hydrody- 
namical details of bubble dynamics. The various dynamical features have been studied 
carefully in 0, where it was argued that field strengths of the order of 10~^^ G on a 
10 Mpc comoving scale could be achieved in EW phase transition by this mechanism 
(and in the case of QCD, even larger fields). As discussed in [§], after the phase tran- 
sition hydromagnetic turbulence is likely to enhance the seed field by several orders of 
magnitude, thus making the primordial field a plausible candidate for the seed field. 

In the present paper we will focus on magnetic fields created in bubble collisions, 
following the treatment of Kibble and Vilenkin who showed that in a collision 
of two bubbles a ring-like magnetic field is formed. (Bubble collisions have also been 
treated in |^ , but mainly with an eye on the defect formation) . The starting point is an 
abelian Higgs model, the properties of which are likely to refiect the properties of the 
full EW SU{2) X U{1) model. In Section 2 we discuss the collision of the bubbles, and 
in particular the multiple collisions, and show that in all cases the resulting magnetic 
field looks qualitatively the same. In Section 3 we introduce diffusion and consider very 
slow collision velocities, which are typical to the EW case. We find out the average 



field by folding in the spectrum of separation of nucleation centers and performing the 
volume average over the randomly inclined ring-like magnetic field configurations. In 
Section 4 we present our conclusions. 

2 Bubble collisions 
2.1 Basic features 

Let us begin by recapitulating some of the features of colliding bubbles in abelian 
Higgs model We shall begin by first ignoring diffusion. The starting point is the 
U(l)-symmetric lagrangian 



A gauge invariant phase difference can be defined in terms of an integral over the 
gradient D^Q [Q. Before the collision, the phase angle within each bubble may be 
taken constant. Following Kibble and Vilenkin 0] we assume that inside the bubble 
the radial mode p settles rapidly to its equilibrium value rj and can thus be treated as a 
constant. It then follows that a Klein-Gordon equation holds for A and separately: 



where X = oi Q. These are not independent but are related by virtue of the 
Maxwell equations: 



where = —ep^{d^Q + eA^). The simplest case is that two bubbles nucleate, one at 
(x, I/, t) = (0, 0, Zq, 0) and the other at (x, y, z, t) = (0, 0, —zq, to), and keep expanding 
with velocity v even after colliding. Their radii at the collision are denoted by Ri and 
J?2- It is easy to find out the intersection volume of the bubbles: 




(1) 




{d^d^- + eW)X = , 



(3) 



d,3' = 0, 



(4) 



x^ + y^ = Rl{t) -{z + zaf ■ z>0 ,x'^ + y^ = Rl{t) - {z - z^f ; z < . (5) 



Denoting Ri{t) = vt and R2{t) = v{t — to) one may solve for the time of intersection: 







In this Section we will always take v = \ and to = for simplicity, so that R\{ti) = 
R2{ti) = R. As we will discuss later, this is not true for EW phase transition where 
V -C 1, but nevertheless the assumption serves a useful illustrative purpose. 

Because of the symmetry in the bubble collision, one can now assume |0] that the 
phase angle G is actually a function of z and r = y/t'^ — x'^ — y'^. Also, the symmetry 
of the problem dictates that in this case the electromagnetic potential has the form 
j^fi — ^fij-^^^ J where / is a function to be determined later. One then finds that in the 
intersection volume the electric and magnetic fields read (in cylindrical coordinates) 

2tr d d 
E = - — g^f{r,z)er-t—f{T,z)e,, 

B = r^fiT,z)e^, (7) 

where r = a/x^ + y'^. 

The solutions for Q and / , in the gauge = and with p =constant, are obtained 
with the initial conditions @] 

Q\r=R = 0oe(^ - Zi) , drQ\r=R = 0, 
f\r=R = , drf\r=R = e{z - Zl) , (8) 



and read 



^ 60-R r°° dk . ^ . , r , ^, 1 . , ^.^ 

B = / —sm.k[z — Zi)[co^uj[T — R)^ — smcu^r — itjj , 

vrr J-00 k uR 

I = ^/ —smk{z - zi)[ cosuj{t-R) 

nr-^ J-00 k oj^R 

+ (- + 4^)sinu;(r-i?)], (9) 

where oo = y/W+~e^ and R is the radius of the bubbles at the collision and zi the 
point of first collision on the z-axis. 

In Fig. 1 and 2 we display the time evolution of the absolute value of the magnetic 
field and O. One can see that B spreads out in the whole intersection region (which 
in Fig. 1 and 2 correspond to the range 50 — t < 2; < 50 + 1) and oscillates rapidly 
with increasing frequency as one moves from the center of the regime to the edge. As 
time goes on, the amplitude of B oscillation decreases as while the frequency 
increases. Therefore the mean field at sufficiently large scales, of the order of 10 /[er]), 
is zero everywhere else except in the middle of the collision regime. The energy density 
1^00 B'^dV/V can be seen to scale as 1/t. 




Figure 1: The time evolution of the phase angle B along the radius r = 1 for t=10, 
20, 30 and 40 after the initial collision. The radius of the bubbles at the collision has 
been chosen here R = 10, and the collision point on the z-axis is zi = 50. (The units 
are such that erj = 1). 

2.2 Full-empty collisions 

After the bubbles have intersected, they will be subject to new collisions. Indeed, first 
order phase transitions are locally completed by merging of several bubbles. Let us 
assume that a sufficiently long time has elapsed since the first collision so that we 
may approximate the result of the collision by a spherical bubble with radius Ri. It 
contains a magnetic field and a phase angle (we shall call it "full" bubble), which are 
given by Eq. (^). Strictly speaking, these solutions cannot be valid inside the whole 
full bubble, but we assume that they are valid in the collision region where the full 
bubble merges with another bubble. 

When a full bubble collides with a recently nucleated bubble with no gauge field 
inside (which we shall call an "empty" bubble), some of symmetries of the empty-empty 
collision remain. Because of the reduced symmetry, we will now take r = \Jt'^ — x^. 

The electromagnetic potential in the full-empty case can be obtained by approxi- 
mating the initial conditions by ,^=o,z>o = /o(l — e~'"^)(l — e'''~^^')'^H{z2 — z), 
where H is the Heaviside theta-function, /q the asymptotic form of the electromagnetic 
potential and drf\T=Ri,y=o,z>o = 0. Z2 is the point of collision on the z-axis, which is 
chosen as the symmetry axis of the collision. These initial conditions should approxi- 
mate quite well the average behaviour of the gauge potential in the 'full' bubble, which 




^ 0.000 --- 




0.0 20.0 40.0 



0.0 100.0 



0.0 20.0 40.0 60.0 



0.0 100.0 





0.0 20.0 40.0 60.0 



0.0 100.0 



0.0 20.0 40.0 60.0 



0.0 100.0 



Figure 2: The time evolution of B along the radius r = 1 for t=10, 20, 30 and 40 
after the initial collision. The radius of the bubbles at the collision has been chosen 
here R = 10, and the collision point on the z-axis is zi = 50. (The units are such that 
ei] = 1). 



subsequently collides with an empty one assuming that the bubbles are big enough 
so that the vector potential A can be taken to point to a constant direction in the 
collision region of the full bubble. 

The vector potential is zero at 2; = and Ai = Xifo in the asymptotic 2; ^ 00 
region, where /o is given by 



/o = f{z 00) 



BovRiR-T ( 1 

cos(e?7(r — R)) + r + 



erjR 



erjR 



sin(er/(r - i?)) . (10) 



The parameter C can be found by matching the z-derivatives of /o(l — e '"^)(1 — 
g-C22)-i ^ ^ ^ 0. 

This approximates rather well the overall behaviour of /, but does not take into 
account the rapid oscillations involved. However, when diffusion is included (Sect. 3), 
it is precisely these rapid oscillations which get diffused. An expression for A, obeying 
the approximate initial conditions, can then be found, but it cannot be expressed in a 
simple form and we do not reproduce it here. A more tractable picture can be obtained 
by making the assumption that the empty bubble is much smaller than the full bubble. 
This also seems a natural assumption. We may then write / = /o + 5/, where 5f is 
a perturbation and /o is the asymptotic solution obtained from the collision of two 
empty bubbles, Eq. (|T0|). 



Let us choose a new coordinate system so that the electromagnetic potential in the 
full bubble can be written in the form 



A = i?i/o sin ae 



y ' 



'111 



where /o is the asymptotic form of /, Eq. (|TOD, a the intersection angle and Ri the 
radius of the large bubble. Also, we take B to be at its asymptotic value as given by 



Or 



-^(cos(e,(r-i?)) + 



(12) 



Here Goo is the original phase difference of the two initial colliding empty bubbles, and 
R refers to the radii of the these bubbles. The initial conditions can now be written 
in the form 



Q\r=Ri = Q0^{Z - Zl) + Ql , dr<S>\r=B., 



0, 



f\r=Rr = a{l - e{z - Zl)) , drf\r=R^ = b{l - e{z - Zl)), (13) 

where a = /o/2 and b = —erf{Qo + ?>fQ/2erf)/Ri. Here r is once again = t^ — x^—y'^ 
and tu^ = fc^ + e^rf. One then finds that 

QqRi 



e 



+ 



TXT 
1 
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cos/(z - Zl) + '^^^^^][cos {JeW^-{r - Ri)) 
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sin( JeV + ^^(^ - ^i))]^^^ 



(14) 



and 

/ 
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au) Ri — cos{2luRi] 
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(15) 



From the complicated expression Eq. (^) one then finds that approximately 

B~i?isin«(9,/)e,. . (16) 

To take into account the angle a between B in the full bubble and the line of collision, 
one should substitute y ycosa — 2; sin a and z 2; cos a + ysina. Qualitatively, 
the situation is however equivalent to the case of empty-empty collision. As can be 
seen from Eq. ([T6|), the resulting magnetic field is again a ring, although it is slanted 
by the angle a with respect to the z-axis. Its time evolution is also similar to what is 
presented in Fig. 1. Because almost every collision that takes place between empty 
and full bubbles can be treated in the approximation above, the structure of magnetic 
fields generated will all be essentially equal. 



2.3 Subsequent collisions 

As time goes on, more and more bubbles may have collided at least once, and the 
subsequent collisions take place between already "full" bubbles, i.e. bubbles which 
contain magnetic fields. In such collisions there are no symmetries left, and the result- 
ing equations of motion are much more complicated compared with the previous cases. 
However, in the case where one of the bubbles is much larger than the other, it is again 
possible to approximate the small bubble as a small perturbation on a bubble with 
infinite radius. We consider a collision where the magnetic fields inside the bubbles 
point to arbitrary directions. The bubbles may again be taken to collide along the 
2-axis with the gauge potentials almost constant in the colliding region. Because all 
symmetries are lost, one can not assume that Ai = Xif. In general, however, these 
collisions can be very well approximated with the initial conditions given by Eq. (|1^) 
for B and A. Therefore, qualitatively the full-full collision does not differ from the 
empty-empty collision discussed in Sect. 2.1. 

Because both B and A obey the Klein-Gordon equation, the solutions are simple 
waves. That is why in the case of two colliding bubbles the qualitative behaviour of the 
generated magnetic field is found not to depend on whether the bubbles had already 
collided or not. A more involved situation arises when one considers a simultaneous 
collision of more than two bubbles. The waves can then interfere and produce more 
complicated patterns. 

To demonstrate this, let us assume that 3 (empty) bubbles move along the 2;-axis 
and collide simultaneously. We can then use symmetry and again define A = x/. The 
initial conditions read as in Sect. 2.1, but now for the three bubbles. We have solved 
the equation of motion numerically, and the results are presented in Fig. 3 and Fig. 
4, where we show the time evolution of B and B. 

A simple way to discuss the collision of several bubble is to assume that one of the 
bubbles is much larger than the others, as was done in Sect. 2.2. When r oo, the 
equation of motion is approximately the wave equation, and thus the small bubbles 
only create small wave-like perturbations which subsequently interfere with each other. 
Our conclusion is that, qualitatively, all possible collisions between the bubbles look 
like collisions of empty bubbles. 

3 Colliding electroweak bubbles 
3.1 Flux spreading and diffusion 

The magnetic field generated in bubble collisions will be imprinted on the background 
plasma. In the early universe electrical conductivity is high but not infinite. When 





Figure 3: The time evolution of in symmetric collision of three bubbles each having 
radius i? = 10 at the collision. Here r = 1 and time refers to the time elapsed after 
the initial collision of the bubbles. The points of initial collision of the bubbles on the 
z-axis are 40 and 60 (and units are 67] = 1). 



1 GeV ^ T < mw it is given by |T0 



a ^ 6.7T 



(17) 



When T ^ mw, and above the electroweak phase transition, one should also account 
for the W, Z and the higgs, but their presence will change a only slightly (note 
that quarks, including the top, contribute very little ||10[ to a). Therefore we will 
adopt Eq. ( p!7| ) as appropriate for electrical conductivity both in the broken and in 
the unbroken phase. Finite conductivity gives rise to a diffusion term in the MHD 
equation 



(5B 

- = Vx(vxB) 



-V X (V X B) 



a 



so that diffusion begins with small scales, and diffusion time td ~ aL^. 

In real cosmological phase transitions the velocity of the bubble wall is definitely 
less than c. In such a situation the magnetic flux will escape the region of intersection 
and penetrate the interior of the bubbles, and eventually the false vacuum outside @]. 
This issue can be addressed by taking the bubble radius to be R(t) = vt and assuming 
r = + ^ 1 so that we may still take r = \/vH'^ — r^. (Here we assume that 
the nucleated bubbles have the same radius.) To include dissipation in the equation 






Figure 4: The time evolution of S in a symmetric collision of three bubbles each having 
radius R — 10 the collision. Here r — 1 and time refers to the time elapsed after 
the initial collision of the bubbles. The points of initial collision of the bubbles on the 
z-axis are 40 and 60 (and units are erj — 1). 




Figure 5: B with = 1, r = 1, cr = 7, i? = 10 and t = 10, 20, 30 and 40 after the 
initial collision. The point of initial collision on the z-axis is zi = 50 (and units are 
er) = 1). 





Figure 6: B with r = 1, i? = 10, a = 7, t = 30 after the initial coUision and v = 0.8, 
0.1, 0.01 and 0.001. The point of initial collision on the z-axis is Zi = 50 (and units 
are erj = 1). 

of motion properly, we add a conduction current term to the Maxwell's equations 
[|] and obtain the following equation of motion for /(r, z): 



]|| + (^-lp(l + .^)-2(.^ + 0} 



r T' 



+ a 



+ e T] J = 0. 



2VT 



dr dz"^ 



(19) 




At the time of intersection, t = R, with R ^ l/T, so that l/crr ~ 1/aR ^ 1. 
Neglecting higher order terms results in 



(20) 



The consequence of the diffusion term in Eq. (|20|) is to smooth out the rapid oscillations 
of B, as is demonstrated in Fig. 5, where we display the numerical solution of Eq. (^) 
for the case v = 1. Fig. 5 should be compared with Fig. 2, which does not include 
diffusion. 

As will be discussed below, in electroweak phase transition the bubbles will in fact 
intersect with non-relativistic velocities. In Fig. 6 we show the velocity effect on the 
created magnetic field. As can be seen from Fig. 6 (and comparing with Fig. 5), the 
smaller the velocity, the smaller the magnetic field. 

Fig. 6 does not yet take fully into account the MHD equation, which should hold 
also outside the intersection region. Taking Eq. (|l^) in conjunction with Eq. (|l^) 
results in the escape of the magnetic flux, which is most easily seen in the case of low 
conductivity (a = 0) and is demonstrated in Fig. 7. In Fig. 8 magnetic diffusion 
is again switched on. Again the magnetic field escapes from the intersection region 
and moves outwards with the speed of light. Here our conclusion is different from 
P], where it was argued that with bubble wall velocities in the range f ~ 0.1 — 1 no 




Figure 7: B with r = 1, = 10, cr = 0, t = 30 after collision and v = 0.8, 0.1, 0.01 
and 0.001. The point of initial collision on the z-axis is zi = 50, and the outer edge 
of the intersection region is at ^ 50 ± vt. The outer edge of the magnetic field is at 
^ 50 ± t (and units are erj = 1). 

magnetic field should be present outside the collision region. For very low velocities 
{v < 0.01) with a = 7 the magnetic field looks however just like the case with a = 0. 
This is due to the fact that for small bubble wall velocity the magnetic field outside 
the bubbles has enough time to decay. Note, however, that decay time is dependent on 
the scale of magnetic structure, i.e. the radii of the bubbles, which in these examples 
are very small compared to the realistic EW case ( the reason being that bubbles with 
very large radii are not easily tractable by numerical analysis and that in the region 
of interest, r ~ i?, the structure of the magnetic field does not depend essentially on 
R). 



3.2 Electroweak phase transition 

The real-time history of a first order electroweak phase transition depends in an essen- 
tial way on the hydrodynamical, non-equilibrium dynamics. Recent lattice simulations 



ri[] have now provided reliably the parameters pertaining the phase transition in the 



case ruH < mw, and a real-time simulation of the bubble growth and collision has 
been performed in [Q. In the following we shall make use of these results. 

When bubbles of broken symmetry are first nucleated, they usually grow rapidly. 
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Figure 8: B along r = 1 with R = 10, a = 7 and t = 30 after collision with v = 0.8 
and 0.1. For the cases v = 0.01 and 0.001 the field looks like in Fig. 6. The point of 
initial collision on the z-axis is zi = 50, and the outer edge of the intersection region 
is at ^ 50 ± vt. The outer edge of the magnetic field is at ^ 50 ± t (and units are 
er] = 1). 

In the electroweak case the initial growth of the bubble wall is by subsonic defiagration, 
with velocities of the order of 0.05c, depending on the assumed friction strength. The 
wall is preceded by a shock front, which may collide with other bubbles. This results 
in reheating, oscillations of the bubble radii, but eventually a phase equilibrium is 
attained. The ensuing bubble growth is very slow and takes place because of the 
expansion of the universe 0. Note that because the universe has been reheated back 
to Tc, no new bubbles are formed during the slow growth phase. 

As an example, let us choose as the reference point the set of numbers related to 
the case of a small friction coefficient 

rave = 9.5 X lO^hn ; v = 1.2 x 10"^ ; 5 = 0.03 , (21) 



where Vave is the average distance of the nucleation centers, v is the velocity of the 
bubbles after the collision of the shock fronts, and 6 is the fraction of the volume in 
the broken phase at the onset of the slow bubble expansion (the initial velocity of the 
wall in this case is Vmit = 0.089). The numbers refer to the case of mn = 68 GeV, 
a weak first order transition. To be definite, we will adopt = 100 GeV and = 
3.55 X 10^3 GeV-\ 

In EW phase transition the collision of the bubbles thus takes place with a very 
slow velocity, and the bubbles are very large. Therefore the discussion in the previous 
Section is not directly applicable, but we may nevertheless assume that the gross 
features hold, in particular that the collision between EW bubbles will (for practical 
purposes) always resemble the collision of two empty bubbles. 

In Fig. 6 we displayed the numerical solution to the equation of motion Eq. (|20|). 
For large EW bubbles numerical methods are not accurate enough to provide a reliable 
result. Therefore we need an analytic esimate of Eq. (^). Let us write f{T,z) = 



6(r)c(z) and use the initial conditions f{R,z) = and drf{R,z) = erfQo ■ e{z)/R, 
where e{z) is the sign-function. Then in the small r = a/xM-^ limit we find that 

ttR 



Sinh(i^/(T2(1 + t;2)2 _ lQq2yi(^j- _ R)/2v'^^ 
/(72(l+t;2)2_ 16^2^4 



cos(fc(2 - 2i)) ^^-^^ ' ' ' ' — -dk , (22) 



where zi is the point of initial collision on the z-axis, = {k'^+e^iff) and = v^t^—r'^. 

From Eq. ([2^ ) one sees that i? at extremely small scales {L < Av/a = 7.1 x 10^^/T) 
decays essentially instantaneously, although strictly speaking one cannot apply MHD 
at length scales less than the interparticle distance ~ 1/T. At scales L > 1/T Eq. (^) 
implies that at a given time > R/v the magnetic field at length scales 

1/2 



has already decayed. Note that at the outer edges of the intersection region diffusion 
cannot remove B, unless r ~ 0, L ^ R. 

Roughly speaking, at the onset of the slow bubble expansion, the bubbles almost 
touch. Once they have collided, the phase transition will be completed in a time 
At ~ 0{R/v). Although strictly speaking Eq. ( p2D is valid only for r <^ vt, let us use 
it to estimate the magnitude of B at the end of phase transition. From Eq. (^) we can 
see that the magnetic field has decayed from the centre of the collision region. Eq. (^) 
tells us that new magnetic field is created only in the region where r ~ i?. The fiux 
escapes, as discussed in Sect. 3.1, and dissipation outside the bubbles is very slow (at 
scales much larger than L in Eq. (^)). Inside the bubbles there is also a region where 
B has not yet decayed. The largest field is obtained, however, around r ~ i?, so that 
the field looks like a narrow ring around the z-axis at z = and r = ^vH'^ — R?. From 
Eq. (|22|) we then obtain an estimate for the strength of the magnetic field in the ring 



i?(i?)^-^^72 + 27i?, (24) 

where we have defined 7 = f At = 0(\)R- This yields the size and the coherence 
length Lo of the field as 



B ~ 4.0^72 + 2-iRlR GeV" = 2.0 x W^^i'^ + l-^RjR G, 
Lo ^ 6.5xlO^((l + ^^^±^)V2_i)i/2/T^, (25) 



where we have assumed Bq = 1, Tc = er/ = 100 GeV and used Eq. 

Thus we may argue that at the end of EW phase transition there are about 
{tu/raveY ~ 10^^ magnetic rings within each horizon volume. To obtain a reliable 



estimate of the average B, we should average over the possible bubble sizes, on which 
B depends. The spectrum of separation of the adjacent shocked spherical bubbles in 
a first order phase transition has been estimated in and reads 

= ^ (-X R'[ewi-S'R/v.)i^ - |) +exp(-2Sm)(f^ + |)] , (26) 
185 \VsJ 2vs 3 4:Vs 3 J 

where Vs = 1/v^ is the velocity of the shock front, and 

S'/vs = S'{tf)/vs = (vr)i/Vr,„, = 4.35 x 10'^ GeV (27) 



is the derivative of the tunneling action at nucleation time t/ |T3|, and the number is 
for the reference values Eq. (|21|). Assuming that Eq. (^) also gives the distribution 
of the bubble radii at collision, we thus arrive at the average magnetic field 



{B) = J P{R)B{R)dR = 3.1 X 10^° G, (2^ 



where we have taken 7 = 2R for definitess. 

Thus, on the average, each volume of a radius Vave contains a large ring-like field 
{B), but with the planes of inclination randomly distributed. For large volumes, we 
should average over all possible inclinations, which corresponds to a random walk on 
the 2d surface of a sphere. This results in a highly entangled field with a root-mean- 
square value (Brms) = \/j2n,k=i{^n) (Bfc)/A^^ = {B)Lo/L. Therefore, the rms-field at 
a given time and comoving scale L reads 

where the factor frit) accounts for the turbulent enhancement of B at large length 
scales, which persists until the plasma becomes matter dominated and even in the 
presence of large plasma viscosity (i.e. Silk damping) The reason for such an 
inverse energy cascade is the non-linear nature of the MHD equations. Effectively, the 
small magnetic loops will merge to form larger magnetic loops, thereby transferring 
energy to larger length scales. Numerical simulations, using so-called shell models to 
simulate the full 3d MHD equations, together with scaling arguments, suggest that 
frit) ~ with p ^ 0.25. The total enhancement may thus be estimated as /j- — 10^ 
(but with somewhat large uncertainty). At the comoving scale of 10 Mpc today we 
thus find for the set of the reference values Eq. (^) (which implicitly assume that 
niH = 68 GeV) 

Brra, ^ lO^^l^ . (30) 

(Here we assumed that the transition from radiation dominated to matter dominated 
era takes place at t ~ 10^ yrs). It is encouraging that the magnetic field given by 
Eq. (|30| ) appears to be of the correct order of magnitude to provide the seed field for 
the galactic dynamo. 



4 Discussion 



Although we have completely neglected the non-abelian nature of the electroweak 
theory, it is unlikely to affect the gross features of true electroweak bubble collisions. 
They, like in the abelian Higgs model we have considered, are bound to depend mostly 
on the existence of a phase difference and magnetic diffusion (although the abelian 
anomaly could be important for magnetic fields [0). We have not treated the full (and 
complicated) hydrodynamics of the bubbles either. These may include deformations 
of the spherical bubbles and kinematic viscosity playing an important part in the 
dynamical evolution. If the surface tension of the electroweak bubbles is not extremely 
low, sphericality should be a good approximation even during bubble collisions. The 
resulting magnetic fields, although not necessarily perfectly ring-like, should be of the 
same order of magnitude as discussed here. For very low surface tension the situation 
might be different. 

Our treatment is obviously a simplification of the highly complex chain of bubble 
collisions. However, it serves to emphasize the decisive role the velocity of the colliding 
bubbles play in magnetic field generation. Velocity itself depends on hydrodynami- 
cal details such as the strength of the shock front, and in the final analysis on the 
parameters of the Higgs potential, in particular on the mass of the higgs. Increasing 
the higgs mass decreases the velocity of the colliding bubbles so that the magnetic 
field would also decrease according to Eq. (|2^). For example, moving from rriH = 50 
GeV to rriH = 68 GeV entails a reduction of f by a factor of ten 0|. At the same 
time the surface tension decreases by a factor of twenty. It is also known []TB[ that for 
high enough higgs mass, the transition is no longer of first order. It therefore seems 
likely that for large enough higgs mass the present treatment is no longer valid, and 
no magnetic field is created. 

The fact that reasonable assumptions seem to produce primordial magnetic fields 
which could serve as the seed field for the galactic dynamo nevertheless emphasizes 
the significance of the electroweak phase transition. There now exist proposals to 
measure the strength of the intergalactic magnetic field to a very high precision [ p^ . 
It is interesting that such measurements could provide information also on the nature 
of the electroweak phase transition. 
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